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Abstract 

We show that large (> 100 GeV) mass splittings between the charged Higgs boson 
(H^) and the neutral Higgs bosons (H° and ^4°) are possible in the Minimal Super- 
symmetric Standard Model (MSSM). Such splittings occur when the a parameter 
is considerably larger than the common SUSY scale, Msusy> an d have significant 
consequences for MSSM Higgs searches at future colliders. 
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1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) |[| is currently the leading can- 
didate for physics beyond the Standard Model (SM). At the tree-level the Higgs sector of 
the MSSM takes the form of a two Higgs doublet model (2HDM), where the coefficients of 
the quartic scalar terms are functions of the electroweak gauge couplings [0 . Five physical 
Higgs bosons are predicted: a charged pair (H + ,H~), two CP-even scalars (h°,H°), and 
a CP-odd A . In the tree-level approximation the lightest Higgs h° satisfies M^o < Mz, 
while Mh± ~ Mho m M40 for Mh± > 200 GeV. Only two parameters are needed to fully 
parametrize the tree-level potential, and these are usually taken to be tan/3 (= 1)2/ v±, 
where Vi is the vacuum expectation value of a Higgs doublet) and one of M H ±, M A o. 

At the 1-loop level the coefficients in the Higgs potential receive contributions from 
virtual particle loops and thus become complicated functions of several (SUSY) param- 
eters. In addition, new quartic scalar terms are generated. This changes the tree-level 
mass relationships, most notably weakening the above mass bound to M h o < 130 GeV . 
The near degeneracy relationship, M#± ~ M^q « Mho for the large (> 200 GeV) mass 
region, is only slightly affected, resulting in mass splittings of (9(10) GeV. Recently there 
has been much interest in the 1-loop effective potential of the MSSM with unconstrained 
CP violating phases [Q, ^j, |6|, [7). It has been shown that the mass splitting between the 
two heavier neutral scalars (now mixed states of CP) can be increased up to 30 GeV for 
large phases. 

In this paper we will show that large mass splittings between the charged and neutral 
Higgs bosons, \M H ± — M A o iH °\ <L 100 GeV, are possible in the CP conserving MSSM in 
a previously ignored parameter space. These large splittings occur when \x is larger than 
the supersymmetric (SUSY) scale (defined as the arithmetic mean of the stop masses) by 
a factor of 4 or more. Although sizeable corrections to the tree-level mass sum rules can 
be implicitly found in previous works e.g. ||, we identify the parameter space for the 
largest splittings and discuss the phenomenological consequences. Importantly we show 
that this parameter space is consistent with current experiments. 

Such splittings would break the commonly assumed degeneracy relation M#± ~ 
Mho ~ M A o for Mh± > 200 GeV, which should no longer be taken as a prediction 
of the MSSM. The Higgs mass spectrum of the MSSM may resemble that of a general 
2HDM or other extended (SUSY) Higgs sectors with scalar singlets etc. Knowledge of 
the maximum possible mass splittings among the Higgs bosons in the MSSM may be 
crucial in distinguishing different models at future colliders. In particular, measurements 
of the mass splittings provide important information on the structure of the underlying 
Higgs sector, especially since Higgs branching ratios (BR) and cross-sections are often 
very similar in many popular models. Moreover, new decays channels involving H^^H 
and A would be open and may possess large branching ratios, thus affecting proposed 
MSSM Higgs search strategies at future colliders. 

Our work is organized as follows. In Section 2 we outline our approach for evaluating 
the mass splittings. Section 3 presents our numerical results, while section 4 contains our 
conclusions. 
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2 Mass splittings in the MSSM 



The most general CP violating Higgs potential of the MSSM may be described by the 
following effective Lagrangian: 



+A 3 ($I$ i)($S$ 2 ) + A 4 ($ t 1 $ 2 )($^ 1 ) + A 5 ($I$ 2 ) 2 + A£($J$i) 2 
+A 6 ($ t 1 $ 1 )($ t 1 $ 2 ) + Ae($|$i)($^$i) + A 7 ($S$ 2 )($;$ 2 ) + A;($ 2 $ 2 )($ 2 $ 1 )(l) 

At tree-level Ai — > A4 are functions of the SU{2) and £7(1) gauge couplings while A5 = 
A6 — A7 = 0. In the 1-loop effective potential all Aj receive sizeable corrections from the 
enhanced Yukawa couplings of the third generation squarks. In particular A 5 7^ A 6 7^ A7 7^ 
and are complex in the presence of SUSY phases in A t ,A^ and \x. Explicit formulae may 
be found in 

The restricted form of the tree-level potential limits the possible mass splittings be- 
tween Mh±,Mho and Ma- It can be shown that M#± and Ma are related by: 
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(A4/2 - Re(X 5 ))v 2 



(2) 



glv/2 and A5 = 0, and the above relationship reduces to the familiar 



At tree-level A 4 
sum rule 

M 2 H± -M 2 A = M 2 W (3) 

With the current LEP bound of M A > 90 GeV one finds M H ± — Ma < 30 GeV, with 
approximate degeneracy for Ma > 200 GeV. 

At the 1-loop level A4 and A5 may be written as (keeping the dominant terms): 
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where X t b is a function of A t ,Ab, Msusy and \i. From now on we will focus on the case 
of fi, A t} b being real, with mention given to the case of complex phases where appropriate. 
It can be seen from eq. (4) that fi, A t or /x, At, > AMsusy would overcome much of the 
suppression from the small coefficients, permitting A4, A5 = 0(1). From eq. (2) this would 
give rise to large mass differences M#± — Ma- 

One can see from eq. (2) that M H ± > Ma requires A4 > A5 while Ma > M H ± requires 
A5 > A4. We note that both terms in A5 are positive definite and the term proportional 
to h\ will dominate unless tan/5 is large (which enhances /i&). In the expression for A4 
there may be both constructive and destructive interference among the various terms and 
so the SUSY correction to the tree-level value may take either sign. 
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For the CP even Higgs bosons one has a mass matrix Ms which when diagonalized 
gives the mass eigenstates h and H°. Note that the 1-loop corrected M.s involves all 
Xi, i = 1 — > 7, and its explicit form may be found in [[|. Of interest to us is the mass 
splitting between M H ± and M H o which may be given approximately by: 



r2 71 ^r2 2 / ^4 - - 2ReA f 



Recently it has been shown that the SUSY radiative corrections to the effective Yukawa 
couplings can be important (for a review see |L(J). Such corrections comprise loops 
involving gluino-sbottom and chargino-stop. The modified hb has already been shown 
to sizeably affect b — > |H|,[p]2|, the effective H^tb coupling and the Higgs decay 



h° — > 66 |L4|, [0|. In our analysis we shall restrict ourselves to tan/5 < 20, which is the 
region where the above corrections to hb have minor impact. 

Note that such large mass splittings occur more naturally in extended versions of the 



MSSM which include a singlet Higgs field in the superpotential e.g. ||16|| . In such models 
eq. (2) is modified to include a term ~ \v 2 , where A = \ijv s (v s is the vacuum expectation 
value of the singlet Higgs field). Thus A = 0(g w ) (i.e. gauge coupling strength) may be 
attained with fi = 0(100) GeV, giving rise to large mass splittings. 



3 Numerical results 

We now present our numerical results for the mass splittings M#± — M#o and Ma — Mh± . 
We take two representative sets of parameters: set (A) gives large M H ± — M H o, and set 
(B) gives large Ma — M H ± . 
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Here all masses are in GeV; Mq 3 ,M^,Mi refer to third generation squark soft SUSY 
breaking masses; A t and Ab are the analogous trilinear couplings. We take /i and A ti b to 
be real in the numerical analysis. For these values of tan (3 the terms ~ h\ dominate the 
terms ~ hf in the expressions for the Aj. We define Msijsy as 



Msusy 



)/2 



(7) 



where t± and ti refer to the lighter and heavier stop eigenstates respectively. For set (A) 



M SUSY ~ 450 GeV and for set(B) M SUS y 
splitting m? 2 — m| 



280 GeV. For both parameter sets the mass 



comfortably satisfies the RG analysis requirement 
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< 0.5 



(8) 
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Note that the large value of \x taken in set (A) and (B) does not imply an unaccept- 
ably light t\, since the corresponding entry in the stop mass matrix is //cot/3, which is 
comfortably suppressed for the assumed values of tan (5. Larger values of tan (3 together 
with large /i would invariably generate values of m~ bi which are lighter than the current 
experimental bounds. We also require that all \ remain in the perturbative region. 

In Fig. 1 (a) and (b), we plot M#± — M#o and Ma — Mjj± as a function of tan/3 for 
/i=4500, 4000, 3500 GeV (from above) and /i=3300, 2800, 2300 GeV (from above), respec- 
tively. For the other parameters we take the choices in (||). 

In both Figs. 1 and 2, the line in thin typeset signifies that the corresponding mass 
for one of t\, b\ and M^o becomes smaller than the present experimental bound. 

In Fig. 1(a), we can understand the behaviour of the splittings M#± — M^o from 
eqns.(4) and fl5|). In eq. (4), increasing tan/? enhances the term ~ hfh^ which is negative 
in this parameter space, and thus reduces the positive A 4 /2. A 5 is also positive and is not 
so sensitive to tan/3 since the top quark Yukawa term dominates. However, 2\ Q / tan/3 
is negative and its modulus becomes smaller as tan/3 increases. The combined effect is 
as follows: for small tan/3 the negative term in eq. (5) cancels the positive terms and 
the splitting is small. As tan/3 increases the positive terms dominate, giving rise to 
large splittings. Further increases in tan/3 reduce the still positive A4, while A5 remains 
approximately constant, which explains the descent of the curve for larger tan/3. We 
have also checked that the splittings are relatively insensitive to variations in Ah since the 
relevant terms are suppressed by hf. One might expect that the splittings could be further 
increased by changing the sign of A t , which makes Xq positive. However it is not easy 
to satisfy the scalar quark mass bound with this choice in the CP conserving scenario. 
Larger values of \x give rise to larger splittings, which is explained in Fig. 2(a) below. 

In Fig. 1(b) one can see that Ma — M^± increases with tan/3, which is explained by 
the fact that A4 is decreasing (see above) and thus enhances the splitting from eq. (2). 
Large \i provides the largest splittings (see below). 

In Fig. 2 (a) we plot M H ± - M H o as a function of u,/M SU sy for A t = 2000, 1900, 1800 
GeV (from above). For the other parameters we take the choices in In Fig. 2 (a), 
we can see that u,/Msusy > 6 allows splittings M H ± — M H o > 80 GeV. Increasing A t 
enhances the positive contributions to A4 and A5, which in turn enhances Mjj± — M^o 
(from eq. (0)). 

In Fig. 2 (b), we plot M A - M H ± as a function of u,/M SUS y for M H ± = 150, 250, 350 
GeV (from above), since it turns out that the splittings are not very sensitive to A t . As 
expected, we can see the decoupling behaviour, with the splitting decreasing for increasing 
M H ±. For low values of h/Msusy the 1-loop corrections are not so large and so one has 
approximately the tree-level result M H ± > Ma- As \xjMsusY increases, A5 is enhanced 
to large positive values. In addition, A 4 is decreased from its positive tree-level value. 
Hence the mass splitting Ma — M H ± increases and exceeds 80 GeV for (i/Msusy ^ 10. 

If CP violating phases are allowed in [i and A t) b then As^ become complex numbers 
in general. The mass eigenstates H® and H® are now mixed states of CP, and the large 
splittings Mh± — M H o or M H o — Mn± are possible in a wider range of parameter space 

which satisfies the mass bounds on and ti (bi). 
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The magnitude of the mass splittings has important consequences for the phenomenol- 
ogy of the MSSM. Our results show that M H ± need not be close in mass to H° and A , and 
the mass differences |M#± — Mh°,a\ may be as large as the analogous values in a general 
2HDM or in other extended Higgs sectors. The relation M#± pa Mjjo^a f° r Mff± — 200 
GeV, which is assumed in existing phenomenological analyses, may be broken in the re- 
gion of large fi/M SUSY . We stress that the large values of \i considered here do not arise 
in popular SUSY models, such as minimal supergravity or gauge mediated models. As 
discussed in 0, in order to be compatible with electroweak symmetry breaking, such large 
values of fi would require the soft SUSY breaking mass parameters m 2 and m\ to be of 
the order of \\x\ and negative, and considerable fine-tuning is necessary. 

If M H ± - M H o > 80 GeV then the 2 body decays -> H°W ± would be open. This 
would offer a new discovery channel for if ± at the LHC, and one which is expected to offer 
a very promising signature. Ref. [17| presented a signal-background analysis showing that 
the background is small, and any model which allows a large BR(i7 ± — > H°W^) would 
provide a very clear signal in this channel. We shall show below that such large BRs are 
possible in the MSSM, and consequently would aid the search for H ± at the LHC. 

If M A - M H ± > 80 GeV then the 2 body decays A -> H+W* would be open. This 
would have important consequences for the process /i + pT — > H ± W If at a muon collider, 
which receives contributions from A , H° mediated s-channel diagrams, and was shown 
to be a promising production mechanism for H ± in the MSSM with M^± pa M^o ^o |TH]. 
Any splittings with Ma > Mh± would enhance the rate for cx(/i + /i _ — ► H ± W T ) compared 
to that in [J18|| . For — Mn± > 80 GeV one could have resonant H ± production via 
fi + pT — > A —>■ H ± W T , which was shown to allow very large cross-sections (> 1 pb) in 
the case of the general 2HDM Rl9| . These possibilities will be pursued in a future article 



20 

In Fig. 3 we plot branching ratios (BRs) for the decay H ± — > H°W ± as a function 
of tan f3. This two body decay is open if M#± — Mho > 80 GeV, and is proportional to 
sin 2 (/3 - a). Also plotted is BR(H ± — > /i°IU ± ) which ~ cos 2 (/? - a). We take the same 
parameters as Fig. 1 (a). The thick (thin) lines represent BR(iJ ± — > H°(h°)W ± ) for 
\i = 4500,4000,3500 (from above), respectively. Note that we do not consider decays to 
SUSY particles H ± — * tb, x ± X° e ^ c - The decays iif — > t\bi are not open for this choice 
of Mh± and — > x^X can be closed by choosing suitable values for Mi,M 2 etc. It can 
be seen from Fig. 3 that BR> 20% is attainable for intermediate values of tan/?, clearly 
showing the importance of this new decay channel. 

The behaviour of the BRs can be qualitatively understood as follows: for low tan/3 
the splittings are not large enough to open the channel and for large tan/3 the newly 
opened channel H — > H^W^ cannot compete with the fermionic decay modes. Thus 
the maximum inpact of H H°W ± decays is for intermediate values of tan /3, which is 
the most problematic region for the search at the LHC. The current search strategies 
utilize — > rv T decays EI], which is most effective for tan/3 > 15, or H — > tb 



decays which cover the regions tan/3 < 3 and tan/3 > 25 p2[ . We can also see that 
BR(i/ ± — > /i°IU ± ), although suppressed in the decoupling limit, can be significant. The 
sizeable BRs shown in Fig. 3 would offer very good detection prospects in these channels 
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In general, the mass splittings presented here may give rise to a large contribution to 
the precisely measured electroweak parameter p. For mass splittings > 20 GeV we find 
that the Higgs contribution is always positive and violates the Sp constraint at the 2 (3) a 
level for \Mjj± — Mh-o^I > 150(200) GeV. The dominant SUSY particle contribution (t — b 
loops) to 5p has the same sign as the Higgs contribution f23|, thus further reducing the 
maximum allowed value of \Mjj± — M H o^\. Therefore the mass splittings |M#± — M H o^\ 
presented here are consistent with the Sp constraints. 



4 Conclusions 

We have shown that very large mass splittings M H ± — M H o^ are possible in the MSSM. 
Such splittings occur in a previously ignored region of the MSSM parameter space and 
may exceed 100 GeV, thus strongly violating the commonly assumed degeneracy relation 
Mjj± fa Mho « Ma- The largest splittings arise for relatively large p > QMsusy- The 
previously neglected 2 body decays H ± —>■ H^W^ and A — » H + W~ , H~W + may become 
the dominant channels, which would have important consequences for Higgs searches at 
future colliders. If Ma — M H ± > 80 GeV, a muon collider could copiously produce H ± at 
resonance. 
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Figure 1: M H ± — M H o (left panel) and M A — M H ± (right panel) as a function of tan (3 for 
H = 4500, 4000, 3500 GeV (from above, left panel) and n = 3300, 2800, 2300 GeV (from 
above, right panel), respectively. The thin lines violate the experimental mass bounds on 
h°, ti or b\. For other parameters, see the text. 
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Figure 2: M H ±—M H a (left panel) and M A — M H ± (right panel) as a function of fi/M S usY 
for A t = 2000, 1900, 1800 GeV (from above, left panel) and M H ± = 150, 250, 350 GeV 
(from above, right panel), respectively. The thin lines violate the experimental mass 
bounds on h°, t 1 or b±. For other parameters, see the text. 
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Figure 3: BR(if ± -> iJ (/i )W^ ± ) (thick (thin) lines) as a function of tan/? for the same 
parameter choice as Fig. 1 (a). 
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